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Abstract 

We present a fully relativistic and gauge invariant framework for calculating the cross sections 
of dilepton production in nucleon-nucleon (NN) collisions which is based on the meson-exchange 
approximation for the A^A^ scattering amplitudes. Predictions of our model are compared with 
those of other covariant models that have been used earlier to describe this reaction. Our results 
are also compared with those of the semiclassical models which are employed to get the input 
elementary cross sections in the transport model calculations of the dilepton production in nucleus- 
nucleus collisions. It is found that cross sections obtained within the semiclassical and quantum 
mechanical models differ noticeably from each other. 

PACS numbers: 25.75.Dw, 13.30.Ce, 12.40.Yx 



1 



I. INTRODUCTION 



A recurring feature of the dilepton (e"'"e~) spectra measured in nucleus-nucleus {AA) 
collisions has been the enhancement (above known sources) in the invariant mass distribution 
of the cross sections in the region of the vector meson and u ) pole mass. This has 
been the case for experiments performed for bombarding energies ranging from as low as 1 
GeV/nucleon (DLS data [1]), through the SPS energies (40 - 158 GeV/nucleon) Q, fl fl to 
the energies employed by the PHENIX collaboration at RHIC (which correspond to invariant 
mass of 200 GeV/nucleon) [5]. The enhancement seen at the SPS energies are understood 
in terms of the modification of the p meson spectral function in the hadronic medium ^. 

However, the large dileptons yields observed in the DLS experiment (in ^^C + ^^C and 
^°Ca + ^°Ca collisions at 1-2 GeV/nucleon beam energies) in the invariant mass (M) range 
from 0.2 GeV to 0.5 GeV, are yet to be explained satisfactorily , Q, [ll| . Indepen- 



dent transport model calculations have been unable to describe these data fully even after 
including contributions from (i) the decay of p and u) mesons which are produced dire ctly 
from the nucleon-nucleon (NN) and pion-nucleon scattering in the early reaction phase [l2| , 
(ii) the in-medium p spectral functions 13], (iii) the dropping p mass with corresponding 



modification in the resonance properties [7|, (iv) an alternative scenario of the in- medium 
effect - a possible decoherence between the intermediate meson states in the vector resonance 
decay 10] . This led to term this discrepancy as "DLS-puzzle" [3, which persists even 
now. 

In order to resolve this unsatisfactory situation the high-acceptance dielectron spectrome- 
ter (HADES) has been built ll| which allows to study the dilepton production in elementary 
proton-proton (pp), proton-deuteron (pd) as well as in proton- nucleus {pA) and AA collisions 
with much wider acceptance region for beam energies up to 8 GeV/nucleon. Unlike the DLS 
experiment HADES also measures the dilepton yields in the quasi-free proton-neutron {pn) 
scattering. The first set of data has already been published [ill, [l^ by this group on ^^C 
+ ^^C collisions at beam energies of 1.0 and 2.0 GeV/nucleon. The remarkable fact is that 
these data agree well with those of the DLS collaboration. Therefore, there is no longer any 
question against the validity of the DLS data and the previous failures to explain them by 
various transport models have to do with problems in the theoretical calculations. 

On the theory side, in a recent HSD transport model calculation [15] it has been shown 
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that if one uses larger cross sections for elementary pp and pn bremsstrahlung processes the 
observed dilepton yields in the relevant invariant mass region for the ^^C + ^^C collisions at 1- 
2 GeV/ nucleon can be reproduced. The support for the enhanced elementary bremsstrahlung 



cross sections comes from the calculations of these processes presented in Ref. 



16| within a 



model which is similar to that used in Ref. ^7i] . Althought the calculations perfomed within 



the two mode 
those of Ref. 



s use the same input parameters yet the cross sections of |l6l| are larger than 



17( 1 by factors of 2-4. It has been further argued in Ref. [l5| that elementary 



bremsstrahlung cross sections larger than those of Ref. 17| were reported already in Ref. 
isl ] within a similar type of model which employes however, realistic T matrices to describe 
the initial nucleon-nucleon (NN) collisions. 

In this paper, we examine the issue of the dilepton production in the elementary A^A^ 



collisions in order to high 

n " 

various models 16|. 



3,0, 



ight and understand the differences seen in the predictions of 
18l | for the corresponding cross sections. This is important for a 
proper theoretical description of the new data on the dilepton production in elementary pp 
and pn collisions which are likely to be announced soon by the HADES collaboration. This 
will also have vital implications for the predictions of the different transport models for the 
HADES dilepton yields in AA collisions. Indeed, in a very recent UrQMD transport model 
analysis of the HADES C + C data [19], it has been shown that even without using the 
enhanced elementary bremsstrahlung cross sections of Ref. [l^ the observed dilepton yields 
can be described fairly well at the 2.0 GeV/nucleon beam energy in the region of relevant 
M values. Although, this theory still underpredicts the data in this mass region at the 1.0 
GeV/nucleon beam energy. 

The major difference between models of Refs. jiol . [itI . [isl lies in the method of implement- 
ing the gauge invariance of the A^A^ bremsstrahlung amplitudes. To investigate this issue, 
we have recalculated the cross sections for the dilepton production in elementary pp and pn 
reactions within a fully relativistic and gauge invariant model which is similar to that used 



in Refs. 
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17l | [to be referred as full quantum mechanical model (FQM)] except for the fact 



that we have used a pseudoscalar (PS) nucleon-nucleon-pion {NNn) vertex instead of the 
pseudovector (PV) one employed by these authors. The reason is that with the PS coupling, 
the contact term (seagull diagram) is not involved in the total Lagrangian which however, 
still remains gauge invariant for bare point like nucleons. With a PV NNtt coupling one 
needs to introduce a contact term to restore the gauge invariance of the total amplitude and 
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FIG. 1: A representative of Feynman diagrams for emission of dilepton in nucleon-nucleon collision 
as considered in this work, (a) denotes emission after NN collisions, (b) before A^A^ collision and 
(c) during NN collision. The box represents an off-shell nucleon or a A isobar. 



different approaches 20(] for constructing this contribution lead to very different results 21|. 
The use of the PS NNn vertex makes our discussions free from this ambiguity. 

In the computation of the amplitudes, strong form factors are introduced to quench the 
contributions from high momenta and to include effects due to the compositeness of the 
nucleon. However, this leads to the loss of gauge invariance. We show here that form 
factors at various vertices (hadronic as well as electromagnetic) can still be implemented 
in our model without loosing gauge invariance. We make a detailed comparison of the 
ingredients and predictions of FQM with those of the models of Refs. 16|, ll8| and examine 
their validity critically. The predictions of the FQM is also compared with those of the 
semiclassical models which are used to obtain the elementary A^A^ bremsstrahlung and delta 
isobar contributions to the elementary dilepton production reactions in transport models to 
calculate the cross sections for this reaction in AA collisions. 



II. FORMALISM 

A representative of the lowest order Feynman diagrams contributing to the dilepton 
production within our model is shown in Fig. 1. The intermediate nucleon or resonances 
can radiate a virtual photon which decays into a dilepton (Figs, la and lb). There are 
also their exchange counterparts. In addition, there are diagrams of similar types where the 
virtual photon is emitted from the nucleon line on the right side. The internal meson lines can 
also lead to dilepton emission (see Fig Ic). To evaluate various amplitudes, we have used 
the same effective Lagrangians for the nucleon-nucleon-meson, resonance-nucleon-meson. 
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nucleon-nucleon-photon and resonance-nucleon-photon vertices as discussed in Ref. 17| with 
the sole exception that for the NNn vertex we have used the PS couphng 

jC^NNtt = igNN-K^NlbT ■ ^■k'^N, (1) 

instead of the PV one. The parameters of the model were taken to be the same as those of 
Ref. |l7| where details of calculations of various amplitudes are given. The on-shell equiva- 
lence of the pseudoscalar and pseudovector couplings makes the parameters independent of 
the choice for the type of the NNtt vertex. 

In the computation of the amplitudes corresponding to diagrams shown in Fig. 1 the 
problem of loosing the gauge invariance arises while using the electromagnetic form factors 
if the complete vertex function for the half-off- shell photon production is not used 
The full vertex is given by 

T^ip',p)=ej2 Y: A''Ff'%W\W;q')Oi'A^ip) (2) 

i=l s,s'=±l 

where operators O are defined as Of = 7^, O2 = —io'^'^q^/2mN and 03 = — with rriN 
being the nucleon mass. A's are projection operators which are given by 

A-(rt - (3) 

with W = a/p^ and s = ±1. In Eq. (2) Fi are the form factors. In these equations p and p' 
denote the four momenta of the incoming and outgoing nucleons, respectively and q is the 
four momentum of the photon. We can relate Fi and F3 by using Ward-Takahashi identity 
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24| 



(WTI) |22| 



g,r^(p',p) = e^-^^^[S-\p') - S-\p)], (4) 

where S{p) is the nucleon propagator. It should be noted that WTI does not pose any 
constraint on the magnetic form factor F2. 

It is easy to see that Eq. (4) will restore gauge invariance for the pn bremsstrahlung 



for the exchange of an uncharged meson (see, e.g. Ref. 22]). The current conservation 
condition (CCC) qfj.T'^ = is satisfied if the contributions of diagrams 1(a) and 1(b) are 
added together. In case of the exchange of charged mesons however, the sum of these two 
diagrams does not vanish and one has to also include the contributions from the diagram 
1(c) to fulfill the CCC. 



In actual calculations the hadronic vertices also contain strong form factors that depend 
on the four momenta of the exchanged mesons. Therefore, for uncharged mesons, it is 
sufficient to multiply all the vertices by the same form factor as that of the hadronic vertices 
in order to keep the gauge invariance. On the other hand, the four momentum of the meson 
changes in Fig. 1(c) for the case of the charged meson exchange and the corresponding meson 
current no longer satisfies the continuity equation even after being multiplied by the same 



hadronic form factor. One needs to multiply this form factor in diagram 1(c) 25|] by an 



additional factor -^(A^) which is obtained by letting F(A^) rnultiply the mesonic current 



and then solving for the continuity equation. This leads to 26 1 



2 2 2 2 

F = 1 + + -rf^, 5 

A'^-qi A^ - q( 

where A,^ is the cut-off parameter and qi and q2 are the four momentum transfers at the left 
and right vertices of graph 1(c), respectively. This result can be interpreted as the photon 
coupling to the regular pion (first term) and to "heavy" pion at the left and right vertices 
(second and third terms, respectively). This way of gauging the strong form factor makes 
it possible to use a given form factor for the meson and a different one for the nucleon but 
still fulfill the WTI. 

As far as form factors for the electromagnetic vertices of the nucleons are concerned, we 
note that for real photons, the gauge invariance mandates Fi = 1, F2 = kn and F3 = 0, 
where is the magnetic moment of the nucleon. For the actual case there is a considerable 
uncertainty in these form factors. F3 is never accessible by experiments since O3 j^^ = for 



any conserved current. As in previous studies |16l . 
electromagnetic form factors for the nucleon. 
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27| we have chosen not to include 



We now compare our model with those of Refs. [la, ll8|]. In Ref. [l8l], instead of the 
one-boson exchange picture of our model, the nucleon-nucleon interaction is included via a 
T-matrix that is based on the Paris potential. However, the nucleon current is not gauge 
invariant in this model. These authors rectify this problem in an ad-hoc manner which may 
not have a microscopic basis. In the model of Ref. |l^ , a pseudovector NNtt coupling has 
been used. With this coupling gauge invariance is preserved with a contact term {NNiT'y 
vertex) added to the total Lagrangian. However, once hadronic form factors are introduced 



the gauge invariance will be violated and a gauge restoration procedure has to be applied 21 1 . 
In Ref. jl^ the same method as described above, has been used for restoring gauge 



D 10 
'10-' 






0.1 0.2 0.3 


0.4 


" 


0.1 0.2 0.3 0.4 




■ 1 ■ 1 ■ 1 1 1 : 

V 2.09 GeV 1 


10^ 




1 1 1 1 1 1 1 1 : 

2.09 GeV 1 




^^-^ PP 


1 


10° 






r 




1 


lO"* 


r 








^ 1 


10-^ 










•\t 


10-^ 








Bremsstrahlung 










r 


'■^ 


10"* 


r 


Bremsstrahlung N -| 









^ _-5 




' 





0.2 0.4 0.6 


0.8 


^" 


0.2 0.4 0.6 0.8 



M (GeV/c ) 



M (GeV/c ) 



FIG. 2: [color online] The invariant mass distribution of the A^A^ bremsstrahlung contributions to 
the dilepton spectra in pp (left) and pn (right) collisions at the beam energies of 1.04 GeV and 2.09 



GeV. Results obtained within our model are shown by solid lines while those of Refs. [18!] and 
by dotted and dashed lines, respectively. 



invariance in the case of pp collisions. No contact graph due to pseudovector coupling appears 
in this case. For the np case however, a contact term is needed. Different approaches for 
constructing the additional current contribution (contact term) to restore gauge invariance 



n 

lead to different types of form factors [20] which can yield quite different results [21|. The 
usual practice is to choose a prescription which provides best agreement with the data. 
However, no comparison is shown with the experimental data in this paper. 



III. RESULTS AND DISCUSSIONS 



In Fig. 2, we show the invariant mass distribution of the pp and pn bremsstrahlung 
contributions to the dilepton spectra at the beam energies of 1.04 GeV and 2.09 GeV. Also 



shown here are results for this reaction as reported in Refs. |16l . 1181 ]. First of all we remark 
that the cross sections calculated in the present work are very similar to those reported in our 
earlier work [17] - their shapes are unchanged while absolute magnitudes of the former are 
slightly larger than those of the latter (by about 10%). However, the cross sections reported 
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FIG. 3: [color online] The invariant mass distribution of the A isobar contribution to the dilepton 



spectra in pp (upper panel) and pn (lower panel) collisions the 



results of our model are shown by full lines while those of Ref. 



seam energy of 1.04 GeV. The 



16l | by dashed line. 



for the pp case in Ref. ISj (to be referred as dJM) are larger than ours for invariant mass 
(M) < 0.25 GeV at the beam energy of 1.04 GeV while they are almost identical for 2.09 
GeV in this mass region. At both the beam energies, dJM results are smaller than ours for 
M > 0.25 GeV. On the other hand, the cross sections of Ref. fl^ (to be referred at KK) are 
larger than our results everywhere in the region of M < 0.6 GeV. An important point to 
note is that there is no overall multiplicative factor that differentiates the results of various 
models. 

Despite using the same diagrams, input parameters and gauge invariance restoration 
procedure, our pp bremsstrahlung cross sections are lower than those of Ref. as can be 
seen in the upper left panel of Fig. 2. Of course, in Ref [l^ a pseudovector NNir vertex has 
been used as compared to the pseudoscalar one employed in this paper. In this context, it 
is worthwhile to note that for the real photon production, the covariant model calculations 



do not depend on the choice of the NNtc coupling (PS or PV) as is shown in Ref. 23]. In 
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FIG. 4: [color online] The calculated dilepton invariant mass distribution for the pp collision at the 
beam energy of 1.04 GeV in comparison to the DLS data. The contribution of the tt^ Dalitz decay 
is also shown here which is the same as that in Ref. 



) uai a. 
. In- 



case of dileptons, different results can arise for two couplings from the magnetic part of the 
NN'y vertex. In fact, in Ref. 29|] it is shown that pp dilepton bremsstrahlung contributions 
obtained with the PV NNn coupling are actually smaller than those calculated with the 
PS one at the beam energy of 2 GeV. The calculations presented in Ref. [l7| also support 
this to some extent. Because for the pn collisions, dJM results are not available, in the right 
panel of Fig. 2 we compare our results with KK cross sections only. We see that latter are 
larger than those of ours by factors ranging between 2-3 at both the beam energies. 

In Fig. 3, we show a comparison of our results and those of Ref ]3] for the invariant mass 
distribution of the A isobar contribution to the dilepton production in pp and pn collisions 
at the beam energy of 1.04 GeV. We note that here too the KK cross sections are larger 
than ours by factors of ~ 2 at smaller values of M even though the two models have used 
the same ingredients and input parameters for this part and there is no ambiguity related 
to gauge invariance as the resonance vertex is gauge invariant by its very construction. 

In Fig. 4, we compare the total cross sections of the dilepton production in pp collisions 
with the DLS data at the beam energy of 1.04 GeV. The cross sections calculated within 
our model are folded with appropriate experimental filter and final mass resolution. The 
folded KK cross sections have been obtained by assuming that the folding procedure does 
not affect the ratios of the unfolded cross sections in the two cases. In this figure we have 
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also shown cross sections for the vr*^ Dahtz decay (vr^ 7e^e ) which are the same as those 



shown in Ref. 



ITj. It is seen that KK cross sections overestimate the DLS data for M < 



0.3 GeV where statistical errors are smaller. The data have larger error bars for M > 0.3 
GeV. In this context the HADES data on the elementary dilepton production reactions are 
expected to be useful because of their low statistical error. 

In the several transport model calculations of dilepton production in the AA collisions, 



the usual practice is to calculate the nuc 
photon approximatio n ( SPA) model 31 



eon bremsstrahlung contributions within a soft 



32| and the delta contribution within a Dalitz 



decay (DDD) model [33|]. The corresponding cross sections are added up to get the total 
elementary dilepton production cross sections. In the SPA model, the radiation from the 
internal lines [Fig. 1(c)] is neglected and the strong interaction vertex is assumed to be 
on-shell (which is reasonable only for small photon energies). This cross section is corrected 
by including a multiplicative factor which is the ratio of the phase space available to the 
two-nucleon system with or without the emission of a dilepton of invariant mass M 321]. 

It is desirable to check the reliability of the semiclassical calculations by comparing them 
with results of a full quantum mechanical model. In the upper panels of Fig. 5, we compare 
the invariant mass distributions of the pn bremsstrahlung cross sections obtained within the 
SPA and FQM approaches at the beam energies of 1.04 GeV and 2.09 GeV. We note that 
the SPA (with corrected phase space) model results agree surprising well in shape with those 
of the FQM. This is however, dependent on the values of the np total cross section used in 
the SPA calculations which is parameterized as (j"^(s) = aimN/{s — 4m^) + a2 mb (see, 
Ref. |32l]) where ai = 18 GeV.mb and a2 = 10 mb and s is the square of the invariant mass 
in the incident channel. 

In the middle panel of Fig. 5, we compare the DDD and FQM A isobar contributions 
to the dilepton spectra calculated within two models. We also show here the individual 
contributions of the post-emission (dashed lines) and pre-emission (dotted lines) graphs to 
the FQM cross sections. We note that at both the beam energies pre-emission diagrams 
contribute substantially to the total FQM cross sections for M > 0.3 GeV. In the DDD 
model these diagrams are not included. In the lower panel of Fig. 5 we compare the total 
cross sections obtained by adding the SPA and DDD contributions (termed as semiclassical) 
with those of the FQM. It is clear that for the important intermediate M values the two cross 
sections differ from each other noticeably. Therefore, care has to be taken in interpreting the 
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FIG. 5: [color online] Comparisons of (i) the SPA and FQM pn bremsstrahlung contributions 
[SPA results with and without phase space correction are shown by dashed and dotted lines, 
respectively] (upper panel), (ii) the DDD (dashed-dotted lines) and FQM A isobar cross sections 
(middle panel); Also shown here are the FQM post-emission (dashed lines) and pre-emission (dotted 
lines) contributions, (iii) the sum of SPA and DDD contributions (dashed line) and the FQM total 
cross sections (lower panel). The FQM results are shown by solid lines. 



transport model results obtained by using the semiclassical elementary dilepton production 
cross sections. 



IV. SUMMARY AND CONCLUSIONS 

In summary, we have presented a fully covariant and gauge invariant model for the dilep- 
ton production in elementary nucleon-nucleon collisions employing a pseudoscalar coupling 
at the nucleon-nucleon-pion vertex. With this coupling, the calculations do not involve the 
kind of gauge invariance related ambiguities that are present in those done with a pseu- 
dovector NNtt coupling. 

We find that similar to the results of Ref. 



lower than those of Ref. [l^ by factors of 2-4 for dilepton invariant mass values below 0.6 



171 ]. our NN bremsstrahlung cross sections are 
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GeV for both pp and pn collisions at 1.04 GeV as well as 2.09 GeV incident energies. On 
the other hand, the pp bremsstrahlung results of Ref. [l^ where realistic T matrices have 
been used to describe the initial NN interaction, are similar to those of our model at the 
beam energy of 2.09 GeV while for 1.04 GeV incident energy they are larger (smaller) than 
our values for dilepton invariant masses smaller (larger) than 0.25 GeV. We stress, however. 



that the current arising from the nucleonic diagram in the model of Ref. [18| is not gauge 

invariant. An important point to note is that there is no overall multiplicative factor to 

differentiate our cross sections from those of Ref. Ksil. It is expected that the results from 

□ 

the HADES group on the dilepton production in pp and pd collisions [30| would provide a 
fresh ground for differentiating between various models. 

Calculations performed with a pseudovector NNn coupling are more appealing because 
this coupling is consistent with the chiral symmetry requirement of the quantum chromo- 
dynamics (QCD) [3J] and also because it leads to negligible contributions from the negative 
energy states (pair suppression phenomena) [35|]. However, in calculations with PV A^A^vr 
vertex the contact terms resulting from different prescriptions of restoring the gauge invari- 
ance will have to be carefully examined. This work is currently underway. 

We found that the total dilepton production cross sections in the elementary A^A^ colli- 
sions as calculated within the semiclassical models and used in most transport calculations, 
differ noticeably from those predicted by the full quantum mechanical model. Therefore, 
quantum mechanical cross sections should be used as input in the transport model calcula- 
tions of dilepton production in nucleus-nucleus collisions in order to interprete their results 
properly. 
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